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ABSTRACT 

We  present  Very  Long  Baseline  Array  observations  of  the  kinematics  of  the  water  masers  associated  with  OH 
12.8—0.9,  the  fourth  member  of  the  so-called  water-fountain  class  of  sources.  We  find  that  the  masers  occupy  two 
distinct  regions  at  the  ends  of  a  bipolar  jetlike  structure  oriented  north-south,  with  the  blueshifted  masers  located  to  the 
north  and  the  redshifted  masers  to  the  south.  The  masers  are  distributed  along  arclike  structures  12-20  mas  across  that 
are  oriented  perpendicular  to  the  separation  axis  with  an  angular  separation  of  ~  1 1 0  mas  on  the  sky.  Our  multiepoch 
observations  show  the  two  maser  arcs  to  be  expanding  away  from  each  other  along  the  axis  of  separation.  The  rela¬ 
tive  proper  motions  of  the  two  maser  regions  is  2.7  mas  yr~*  (^105  km  s_1  at  the  assumed  distance  of  8  kpc).  The 
measured  radial  velocity  difference  between  the  northern  blueshifted  masers  and  the  southern  redshifted  masers  is 
48.4  km  s_1.  The  radial  velocity,  when  combined  with  the  proper  motion,  yields  a  three-dimensional  expansion 
velocity  of  58  km  s_1  and  an  inclination  angle  of  24°  for  the  jet.  By  combining  our  radial  velocities  with  historical 
values,  we  estimate  the  three-dimensional  acceleration  of  the  masers  to  be  ^0.63  km  s~ 1  yr~ 1  and  a  dynamical  age  for 
the  collimated  outflow  of  ^90  yr. 

Subject  headings:  circumstellar  matter  —  masers  —  stars:  AGB  and  post-AGB  — 
stars:  individual  (OH  12.8—0.9)  —  stars:  mass  loss 


1.  INTRODUCTION 

In  the  study  of  stars  and  stellar  evolution,  certain  stages  in  the 
stellar  life  cycle  remain  a  mystery  simply  because  of  their  transi¬ 
tory  nature.  The  evolutionary  stage  between  the  end  of  the  asymp¬ 
totic  giant  branch  (AGB)  and  the  planetary  nebula  ( PN )  phase  is 
one  such  example.  As  stars  evolve  up  the  AGB,  they  lose  mass  at 
an  increasing  rate  and  spherically  symmetric  circumstellar  enve¬ 
lopes  (CSEs)  are  formed.  Despite  the  spherical  symmetry  of  the 
mass-loss  process,  a  large  fraction  (~75%)  of  planetary  nebulae 
(PNs)  show  aspherical  (i.e.,  elliptical,  bipolar,  quadrupolar)  mor¬ 
phologies  (Manchado  et  al.  2000).  The  shaping  of  asymmetrical 
PNs  must  therefore  occur  over  the  relatively  short  time  period  of 
1 03  - 1 04  yr  that  the  star  spends  in  a  post-AGB  or  proto-planetary 
nebula  (PPN )  phase  (Kwok  1993).  On  the  basis  of  comparative 
studies  of  compact  and  extended  PNs,  Aaquist  &  Kwok  (1996) 
find  that  the  lack  of  difference  in  the  morphologies  between  the 
two  groups  suggests  that  PN  morphologies  are  primarily  inherited 
from  the  AGB  evolutionary  stage.  Sahai  &  Trauger  (1998)  sur¬ 
veyed  a  number  of  young  PNs  with  the  Hubble  Space  Telescope 
Wide  Field  Planetary  Camera  2  and  found  that  the  majority  of 
them  were  characterized  by  multipolar  bubbles  distributed  in  a 
point-symmetric  fashion  about  the  central  star.  They  also  found 
collimated  radial  structures  and  bright  equatorial  structures  in¬ 
dicating  the  presence  of  jets  and  disks/tori  in  some  objects.  These 
observations  led  Sahai  &  Trauger  ( 1 998)  to  propose  a  mechanism 
in  which  high-speed  collimated  outflows  operate  during  the  late 
AGB  or  early  PPN  stages  of  stellar  evolution.  During  this  phase, 
the  jetlike  outflows  carve  out  an  imprint  in  the  spherical  AGB  wind, 
and  it  is  this  imprint  that  provides  the  morphological  signature 
necessary  for  the  development  of  aspherical  PNs. 

There  is  a  small  number  of  evolved  objects  that  have  been 
shown  to  exhibit  the  types  of  collimated  jetlike  outflows  described 
above.  These  objects,  dubbed  “water-fountain”  nebulae,  are 


thought  to  be  stars  entering  the  post-AGB  or  PPN  evolutionary 
stages.  They  exhibit  both  H20  and  OH  maser  emission;  how¬ 
ever,  the  relative  characteristics  of  the  two  maser  species  differ 
from  those  of  the  typical  AGB  star.  These  differences  are  ap¬ 
parent  in  both  the  spectral  profiles  of  the  two  types  of  masers  and 
in  the  unique  spatial  morphologies  of  the  masers  as  measured  by 
radio  interferometry.  Prior  to  our  study  of  OH  12.8—0.9  (Boboltz 
&  Marvel  2005),  there  were  three  confirmed  water-fountain 
sources:  IRAS  16342—3814  (Sahai  et  al.  1999;  Morris  et  al. 
2003),  IRAS  19134+2131  (Imai  etal.2004),  andW43A(Imai 
et  al.  2002). 

Spectrally,  these  water-fountain  sources  exhibit  double-peaked 
OH  and  H20  maser  profiles  with  the  peaks  symmetrically  dis¬ 
tributed  about  the  radial  velocity  of  the  star.  However,  unlike  the 
typical  evolved  star,  the  H20  maser  peaks  have  a  greater  spread 
in  velocity  than  the  corresponding  OH  masers.  For  IRAS  1 6342— 
3814,  IRAS  19134+2131,  and  W43A,  these  velocity  ranges  for 
the  H20  masers  are  259, 132,  and  I  SO  km  s  *,  respectively  ( Likkcl 
et  al.  1992).  In  addition,  the  expansion  velocities  (>60  km  s  1 ) 
implied  by  the  maser  spectra  are  much  greater  than  the  expan¬ 
sion  velocity  for  the  circumstellar  wind  of  a  typical  AGB  star 
(5-30  km  s~'). 

Spatially,  the  H20  masers  associated  with  these  sources  are 
known  to  trace  bipolar  jetlike  outflows;  hence  the  “water  foun¬ 
tain”  name.  In  the  first  very  long  baseline  interferometry  ( VLBI) 
study  of  the  H20  masers  toward  W43  A,  Imai  et  al.  (2002)  showed 
that  the  water  masers  are  formed  in  a  collimated  precessing  jet 
with  an  estimated  three-dimensional  (3D)  outflow  velocity  of 
145  kms-1.  Interferometric  studies  of  IRAS  16342— 3814  (Morris 
et  al.  2003;  Claussen  et  al.  2004)  and  IRAS  19134+2131  (Imai 
et  al.  2004)  also  found  the  H20  masers  to  exhibit  bipolar  distri¬ 
butions.  The  3D  outflow  velocities  for  IRAS  16342—3814  and 
IRAS  19134+2131  are  185  and  130  km  s_1,  respectively.  The  bi¬ 
polar  jets  traced  by  the  H20  masers,  presumably  along  the  polar 
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axis  of  the  star,  may  represent  the  onset  of  the  axisymmetric 
morphologies  that  typify  PNs.  The  dynamical  ages  of  the  jets  for 
IRAS  1 9 134+2 131  and  W43  A  are  estimated  to  be  ~50  and  ~40  yr, 
respectively  (Imai  et  al.  2002,  2004).  That  for  IRAS  16342-3814 
is  estimated  to  be  ^150  yr  (Claussen  et  al.  2004). 

The  enigmatic  source  OH  12.8—0.9  was  first  classified  by 
Baud  et  al.  (1 979)  as  a  type  II OH/IR  star  on  the  basis  of  its  char¬ 
acteristic  double-peaked  1612  MHz  OH  maser  profile.  The  dis¬ 
tance  to  the  source  is  unknown,  although  Baud  et  al.  (1985) 
included  it  in  a  sample  of  OH/IR  stars  associated  with  the  Ga¬ 
lactic  center.  OH  12.8—0.9  has  been  linked  to  the  infrared  source 
IRAS  18139—1816,  which  is  ~26"  away  (te  Lintel  Hekkert  et  al. 
1989).  The  S1MBAD  astronomical  database,  however,  still  lists 
two  separate  entries  for  the  source,  one  for  each  identifier.  Clas¬ 
sification  of  IRAS  1 8 139— 1816  based  on  IRAS  Low  Resolution 
Spectrometer  (8-23  /.tm)  data  was  performed  by  Kwok  et  al. 
(1997),  who  placed  the  source  in  category  I,  a  group  with  noisy 
or  incomplete  spectra.  Kwok  et  al.  ( 1 997)  also  noted  that  IRAS 
18139—1816  is  a  “25  pm  peaker”:  a  source  whose  IRAS  flux 
density  at  25  pm  is  greater  than  its  flux  density  at  both  12  and 
60  pm.  Other  25  pm  peakers  include  young  PNs  and  carbon 
stars  with  circumstellar  silicate  emission  features  ( Kwok  et  al. 
1997).  OH  12.8—0.9  has  been  also  observed  with  the  PHT-S 
(2.5-11.6  pm)  spectrophotometer  on  board  the  Infrared  Space 
Observatory  (ISO).  From  the  spectra,  Hodge  et  al.  (2004)  clas¬ 
sify  the  source  in  the  4/5. S A  category.  Sources  in  the  4/5  group 
peaklongwardofthe  11.6  pm  limit  of  the  PHT-S.  Objects  in  the 
4/5. S A  subgroup  additionally  exhibit  a  deep  10  pm  absorption 
feature  and  absorption  features  from  H20  at  3  and  6  pm,  C02  at 
~4.3  pm,  and  CO  at  ~4.6  pm  (Hodge  et  al.  2004). 

H20  maser  emission  from  OH  12.8—0.9  was  first  detected  by 
Engels  et  al.  (1986),  who  noted  the  fact  that  the  H20  emission 
peaks  are  outside  the  range  of  the  OH  maser  emission.  Gomez 
et  al.  (1994)  used  the  Very  Large  Array  ( VLA)  to  show  thatthe 
“anomalous”  H20  and  OH  maser  emission  was  spatially  co¬ 
incident  to  within  1"  and  therefore  belonged  to  the  same  source. 
Gomez  et  al.  (1994)  also  found  double-peaked  profiles  for  both 
the  OH  and  the  H20,  with  the  peaks  of  the  OH  separated  by 
~23  km  s~ 1  and  the  H20  peaks  separated  by  nearly  twice  this 
amount,  ^42  km  s  1 .  Although  the  velocity  range  for  the  H20  is 
not  as  wide  as  the  other  water-fountain  sources,  Engels  (2002) 
found  that  the  shape  and  variations  in  the  spectra  were  consistent 
with  the  water-fountain  class  and  compatible  with  an  axisym¬ 
metric  wind. 

In  Boboltz  &  Marvel  (2005)  we  reported  on  our  initial  Very 
Long  Baseline  Array  (VLB  A)  observations  of  the  H20  masers 
associated  with  OH  12.8—0.9.  We  found  that  the  H20  masers 
trace  the  bipolar  morphology  typical  of  the  water- fountain  class. 
This  bipolar  structure  combined  with  the  spectral  characteristics 
of  the  H20  and  OH  maser  emission  led  us  to  propose  that  OH 

12.8— 0.9  is  a  fourth  member  of  the  rare  water-fountain  class  of 
objects.  In  this  article,  we  present  additional  multiepoch  observa¬ 
tions  made  with  the  VLBA  with  the  goal  of  measuring  the  ki¬ 
nematics  of  the  H20  masers  that  trace  the  jet  of  OH  12.8—0.9. 

2.  OBSERVATIONS 

We  observed  the  22.2  GHz  H20  maser  emission  from  OH 

12.8- 0.9  ( a  =  18h16m49.23s,  6  =  -18°15'01.8",  J2000.0)  us¬ 
ing  the  10  stations  of  the  VLBA.  The  VLBA  is  operated  by  the 
National  Radio  Astronomy  Observatory  (NRAO).1  VLBA  spec- 

1  The  National  Radio  Astronomy  Observatory  is  a  facility  of  the  National 
Science  Foundation  operated  under  cooperative  agreement  by  Associated  Uni¬ 
versities,  Inc. 


tral  line  observations  occurred  on  2004  June  21,  2005  July  21, 
2005  November  2,  and  2006  February  24.  Each  observing  run 
was  approximately  5  hr  in  length.  A  reference  frequency  of 
22.23508  GHz  was  used  for  the  H20  maser  transition.  Data 
were  recorded  in  dual  circular  polarization  using  two  8  MHz 
(112.6  km  s~*)  bands  centered  on  the  local  standard  of  rest 
(LSR)  velocity  of  —58.0  km  s-1. 

For  an  unknown  reason,  the  observations  made  on  2005  July 
21  did  not  detect  the  H20  masers  from  OH  12.8—0.9.  Following 
this  nondetection,  we  requested  a  short  (5  minute)  observation  of 
OH  12.8—0.9  with  the  Green  Bank  Telescope  (GBT)  in  order  to 
decide  whether  the  masers  would  still  be  detectable  in  the  pend¬ 
ing  two  epochs  of  VLBA  observations.  On  2005  September  13 
the  GBT  observed  the  H20  maser  spectrum  and  found  masers  at 
levels  of  2-3  Jy.  We  therefore  continued  with  remaining  two 
scheduled  VLBA  epochs  and  subsequently  managed  to  detect 
and  image  the  masers.  The  remainder  of  this  article  relates  only 
to  the  three  epochs  in  which  the  maser  observations  were  suc¬ 
cessful,  hereafter  denoted  as  epoch  1  (2004  June  21),  epoch  2 
(2005  November  2),  and  epoch  3  (2006  February  24). 

The  data  were  correlated  at  the  VLBA  correlator  operated 
by  NRAO  in  Socorro,  NM.  Auto-  and  cross-correlation  spectra 
consisting  of  512  channels  with  channel  spacings  of  15.63  kHz 
(~0.22  km  s_1)  were  produced  by  the  correlator.  Calibration  of 
each  of  the  three  epochs  was  performed  in  accordance  with  stan¬ 
dard  VLBA  spectral  line  procedures  using  the  Astronomical 
Image  Processing  System  (A1PS)  maintained  by  NRAO.  The 
calibration  of  epoch  1  is  described  in  Boboltz  &  Marvel  (2005), 
and  the  two  subsequent  epochs  were  calibrated  in  a  similar  man¬ 
ner.  For  each  epoch,  residual  delays  due  to  the  instrumentation 
were  corrected  by  performing  a  fringe  fit  on  the  continuum  cal¬ 
ibrator  (J 175 1+0939)  scans.  Residual  group  delays  for  each  an¬ 
tenna  were  determined  and  applied  to  the  target  data. 

The  bandpass  response  was  determined  from  scans  onJ1751+ 
0939  and  was  used  to  correct  the  OH  12.8—0.9  data.  The  time- 
dependent  gains  of  all  antennas  relative  to  a  reference  antenna 
were  determined  by  fitting  a  total  power  spectrum  (from  the  ref¬ 
erence  antenna  with  the  target  source  at  a  high  elevation)  to  the 
total  power  spectrum  of  each  antenna.  The  absolute  flux  density 
scale  was  established  by  scaling  these  gains  by  the  system  tem¬ 
perature  and  gain  of  the  reference  antenna.  Errors  in  the  gain  and 
pointing  of  the  reference  antenna  and  atmospheric  opacity  varia¬ 
tions  contribute  to  the  error  in  the  absolute  amplitude  calibration, 
which  is  accurate  to  about  1 5%— 20%. 

Residual  fringe  rates  were  obtained  by  fringe-fitting  a  strong 
reference  feature  in  the  spectrum  of  OH  12.8—0.9.  For  each 
epoch  we  used  the  same  strong  feature  at  channel  velocity  FLsr  = 
—81.6  km  s~*.  The  resulting  fringe-rate  solutions  were  applied 
to  all  channels  in  the  spectrum.  An  iterative  self-calibration  and 
imaging  procedure  was  then  performed  to  map  the  emission  in 
the  reference  channel  for  each  epoch.  The  resulting  residual  phase 
and  amplitude  corrections  from  the  reference  channel  were  then 
applied  to  all  channels  in  the  8  MHz  band.  All  of  the  above  cal¬ 
ibrations  were  then  applied  to  the  data  prior  to  imaging. 

Imaging  of  the  epoch  1  data  is  discussed  in  Boboltz  &  Marvel 
(2005),  and  a  similar  methodology  was  employed  for  the  two 
subsequent  epochs.  For  epochs  2  and  3,  full-resolution  images  of 
2048  x  2048  pixels  (~  1 60  x  1 60  mas)  were  generated  using  syn¬ 
thesized  beam  sizes  of  1.08  x  0.39  mas  for  epoch  2  and  0.88  x 
0.33  mas  for  epoch  3.  Images  were  produced  for  all  spectral  chan¬ 
nels  from  —88.3  kms-1  to  —25.1  kms-1,  forming  an  image  cube 
of  30 1  image  planes.  Since  peak  maser  flux  densities  were  on  the 
order  of  1  Jy,  channel  images  are  noise-limited  and  do  not  suf¬ 
fer  from  dynamic  range  limitations,  as  is  often  the  case  in  maser 
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Fig.  1 . — H20  maser  emission  toward  OH  12.8—0.9  from  epoch  2  (2005  November  2).  In  (a),  ( b ),  and  (c),  the  upper  subpanels  show  the  spectra  formed  by  plotting 
maser  component  flux  density  vs.  LSR  velocity,  color-coded  according  to  maser  velocity.  A  dark  solid  line  in  each  upper  subpanel  represents  the  scalar-averaged  cross¬ 
power  spectrum  on  the  Los  Alamos-Pie  Town  VLBA  baseline.  A  dotted  line  in  the  upper  subpanel  of  (a)  represents  the  scalar-averaged  cross-power  spectrum  of  the 
OH  masers  (with  the  flux  density  scaled  by  a  factor  of  0.5)  from  the  Los  Alamos-Pie  Town  baseline.  The  lower  subpanels  in  ( a ),  ( b ),  and  (c)  plot  the  spatial  distribution 
of  the  H20  masers,  with  point  color  representing  the  corresponding  velocity  bin  in  the  spectrum  and  point  size  proportional  to  the  logarithm  of  the  maser  flux  density. 
(a)  All  H20  maser  components  from  our  VLBA  observations.  A  dashed  line  represents  the  axis  of  separation  between  the  centers  of  the  blue-  and  redshifted  masers  at  a 
position  angle  of  2.6°  east  of  north.  The  cross  represents  the  midpoint  of  the  bipolar  distribution.  ( b )  and  ( c )  Expanded  views  of  the  blue-  and  redshifted  maser  features, 
respectively.  Errors  in  the  positions  of  the  features  are  smaller  than  the  data  points  for  all  panels. 


observations.  Typical  1  a  off-source  noise  estimates  for  indi¬ 
vidual  image  planes  were  8-10  mJy  beam-1  for  epoch  2  and 
10-12  mJy  beam-1  for  epoch  3.  The  analysis  and  extraction  of 
relevant  information  from  the  image  cubes  is  described  below. 

3.  RESULTS 

In  order  to  identify  and  extract  maser  component  parameters, 
two-dimensional  Gaussian  functions  were  fitted  to  the  emission 
in  the  image  plane  of  each  spectral  (velocity)  channel  using  the 
AIPS  task  SAD.  Image  quality  was  assessed  using  the  off-source 
rms  noise  in  the  image.  A  cutoff  flux  density  was  set  at  5  times  the 
rms  noise  in  the  image  plane  containing  the  emission.  Emission 
features  with  flux  densities  greater  than  the  cutoff  were  fitted 
with  Gaussians  to  determine  maser  component  parameters.  The 
errors  in  right  ascension  and  declination  of  the  Gaussian  fits  to 
the  identified  components  were  computed  by  the  AIPS  task  SAD 
following  the  method  outlined  by  Condon  (1997).  These  posi¬ 


tion  errors  ranged  from  5  pas  for  features  with  high  signal-to- 
noise  ratios  to  1 1 0  pas  for  features  with  lower  signal-to-noise 
ratios. 

Figures  1  and  2  show  the  spectral  ( upper  subpanels)  and 
spatial  {lower  subpanels)  distributions  of  the  FLO  masers  toward 
OFI  12.8—0.9  from  the  analysis  of  our  VLBA  images  for  epoch  2 
(Fig.  1)  and  epoch  3  (Fig.  2).  A  similar  plot  for  epoch  1  was  pre¬ 
sented  in  Boboltz  &  Marvel  (2005).  In  the  figures,  panel  a  shows 
the  entire  range  of  FLO  maser  emission  from  the  star.  Panels  b 
and  c  show  enlarged  views  of  the  blueshifted  masers  to  the  north 
and  the  redshifted  masers  to  the  south,  respectively. 

In  the  upper  subpanels  in  Figures  1  a  and  2a,  the  H20  maser 
spectra  for  epochs  2  and  3  are  plotted,  along  with  the  spectrum  of 
the  1612  MHz  OH  masers  from  previous  observations  (Boboltz 
&  Marvel  2005)  as  a  reference.  The  velocities  of  the  OH  maser 
peaks  have  not  changed  significantly  over  time  when  compared 
to  previous  observations  by  Gomez  et  al.  (1994).  We  observed 
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Fig.  2. —  Same  as  Fig.  1,  but  for  epoch  3  (2006  February  24). 


the  velocities  of  the  blue-  and  redshifted  OH  peaks  to  be  —68.0 
and  —43.7  km  s  respectively,  with  a  velocity  resolution  of 
0.18  km  s  1 .  By  comparison,  the  velocities  of  the  two  peaks 
reported  in  Gomez  et  al.  (1994)  are  essentially  the  same  at  —67.4 
and  —44.0  km  s,  to  within  the  velocity  resolution  of  their  VLA 
observations  (1.1  km  s  ').  As  in  epoch  1 ,  the  H20  masers  form  a 
double-peaked  profile  with  a  velocity  extent  greater  than  that  of 
the  double-peaked  OH  maser  profile.  On  the  blueshifted  side  of 
the  emission,  the  peak  flux  density  of  epoch  2  remained  the  same 
as  that  for  epoch  1  at  1.3  Jy  beam-1  and  decreased  slightly  in 
epoch  3  to  0.97  Jy  beam-1.  For  the  redshifted  emission,  the  peak 
flux  density  steadily  decreased  from  f  .0  Jy  beam- 1  in  epoch  1  to 
0.6  Jy  beam-1  in  epoch  2  to  0.2  Jy  beam-1  in  epoch  3.  The  ve¬ 
locity  distributions  for  the  masers  in  epochs  2  and  3  are  very 
similar  to  those  in  epoch  1  and  will  be  discussed  in  detail  in  §  3 .2. 

The  lower  subpanels  of  Figures  lu  and  2 ci  show  the  entire 
spatial  extent  of  the  H20  masers  for  epochs  2  and  3,  respectively. 
Enlarged  views  of  the  north  and  south  regions  are  also  shown  in 
panels  b  and  c  of  each  figure.  We  find  that  in  epochs  2  and  3  the 
overall  spatial  structure  of  the  H20  masers  has  remained  essen¬ 
tially  unchanged  from  that  in  epoch  1 .  There  are  two  distinct  re¬ 
gions  oriented  roughly  north-south  on  the  sky  that  are  separated 
by  ~  1 1 0  mas.  The  position  angle  of  the  axis  of  separation  be¬ 


tween  the  mean  position  centers  of  the  northern  and  southern 
maser  regions  is  ~2.6°  east  of  north.  The  masers  in  the  two  re¬ 
gions  are  arranged  in  arclike  structures  oriented  roughly  east- 
west,  with  the  exception  of  the  feature  located  near  (15.0,-10.0) 
in  Figures  lb  and  2b.  This  feature,  although  apparent  in  the 
epoch  1  image  cube,  did  not  meet  all  of  the  component  selection 
criteria  in  our  initial  analysis  of  the  data  (Boboltz  &  Marvel 
2005).  The  feature  is  discussed  further  in  §  3. 1  below.  The  blue- 
and  redshifted  arcs  of  masers  are  approximately  20  and  12  mas 
across,  respectively,  with  corresponding  opening  angles  for  the 
jet  of  1 3  20  . 

3.1.  Maser  Proper  Motions 

In  Figures  1  and  2,  we  plotted  the  maser  emission  identified  in 
every  spectral  channel;  thus,  an  individual  maser  feature  will 
appear  in  multiple  spectral  channels  and  will  consist  of  multiple 
points  in  the  spatial  distributions  shown  in  the  figures.  In  order  to 
track  and  estimate  the  motions  of  the  maser  features,  it  is  neces¬ 
sary  to  determine  average  values  in  right  ascension,  declination, 
and  velocity  for  each  feature.  These  average  values  were  com¬ 
puted  using  a  flux-density -squared  weighted  average  for  emis¬ 
sion  identified  in  two  or  more  channels  and  spatially  coincident 
to  within  1  mas.  The  flux  assigned  to  the  maser  feature  was  simply 
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TABLE  1 

H20  Maser  Characteristics  Derived  from  the  Epoch  1  (2004  June  21)  VLB  A  Data 


i>lsr  V  <Jsv  Relative  R.A.  or.a.  Relative  Decl.  ODeci. 

(kms-1)  (Jy  beam  " 1 )  (Jy  beam-1)  (mas)  (mas)  (mas)  (mas) 


Blueshifted  Masers 

-85.07 .  0.044  0.009  5.869  0.038  2.520  0.063 

-83.32  .  0.052  0.009  14.806  0.034  -10.133  0.049 

-81.71  .  1.297  0.012  0.000  0.006  0.000  0.008 

-81.67 .  0.123  0.012  -1.813  0.031  -0.322  0.045 

-80.88  .  0.261  0.009  7.647  0.009  0.376  0.014 


Redshifted  Masers 

-34.61  .  0.247  0.009  1.116  0.012  -108.821  0.019 

-33.92  .  0.139  0.012  -5.773  0.018  -108.460  0.027 

-33.72  .  0.498  0.012  2.235  0.018  -108.856  0.030 

-33.27 .  1.041  0.012  -2.807  0.005  -108.597  0.007 

-32.58  .  0.499  0.009  5.297  0.010  -106.906  0.017 


the  peak  flux  in  the  channels  spanned  by  the  feature.  Char¬ 
acteristics  of  the  masers  for  the  three  epochs  are  presented  in 
Tables  1-3. 

Since  we  did  not  use  the  technique  of  phase-referencing  for 
our  three  VLB  A  epochs,  the  absolute  position  of  the  phase  center 
in  each  image  cube  is  unknown.  However,  in  the  reduction  and 
analysis  of  each  epoch  of  observations,  we  used  what  we  believe 
to  be  the  same  maser  feature  as  a  phase  reference.  This  feature 
was  the  strongest  feature  in  each  of  the  three  epochs.  In  addition, 
it  maintained  the  same  velocity  and  spatial  location  relative  to 
the  other  three  northern  blueshifted  masers  common  to  all  three 
epochs.  The  averaging  of  the  maser  identifications  over  several 
channels  resulted  in  slight  (<0. 1  mas)  position  shifts  for  the  ref¬ 
erence  feature  in  each  of  the  three  epochs.  The  coordinate  frames 
for  all  three  epochs  were  therefore  realigned  such  that  the  ref¬ 
erence  feature  (the  third  blueshifted  feature  in  Table  1  and  the 
second  blueshifted  feature  in  Tables  2  and  3)  coincided  with  the 
origin. 

Relevant  maser  characteristics  determined  from  the  analysis 
are  plotted  in  Figures  3  and  4.  Figure  3  shows  the  blueshifted 
H20  masers  in  the  north  of  OH  12.8—0.9.  In  this  figure,  the 
alignment  of  the  reference  features  for  the  three  epochs  at  the 
origin  is  apparent.  Also  common  to  all  three  epochs  are  features 
near  (—2.0,  0.5),  (8.0,  0.2),  and  (15.0,  —10.0).  As  mentioned 


earlier,  the  feature  near  (15.0,  —10.0)  did  not  appear  in  Boboltz 
&  Marvel  (2005)  because  it  did  not  meet  all  of  the  selection 
criteria  at  the  time.  However,  with  the  additional  data  from  ep¬ 
ochs  2  and  3  and  some  reprocessing  of  the  epoch  1  image  cube, 
we  were  able  to  extract  the  relevant  parameters  for  this  feature. 
From  Figure  3  it  is  apparent  that  any  position  differences  be¬ 
tween  the  three  epochs  are  small,  thus  indicating  little  motion 
of  the  northern  features  relative  to  the  reference  feature  over  the 
course  of  the  three  epochs.  This,  however,  is  not  the  case  for  the 
southern  maser  features.  Figure  4  shows  the  redshifted  FLO  ma¬ 
sers  in  the  south  of  OH  12.8—0.9.  Here  there  are  five  features  that 
are  common  to  all  three  epochs.  Clearly  there  is  a  consistent 
proper  motion  of  all  of  the  features  of  nearly  due  south  relative  to 
the  blueshifted  reference  feature.  The  average  angular  separa¬ 
tions  of  the  north  and  south  masers  are  107,  110,  and  111  mas  for 
epochs  1, 2,  and  3,  respectively.  As  discussed  in  Boboltz  &  Marvel 
(2005),  the  distance  to  OH  12.8—0.9  is  not  well  known.  Best  es- 
timatesplace  the  object  near  the  Galactic  center  at  £>  ss  8kpc(Baud 
et  al.  1985).  If  we  assume  this  distance,  then  the  linear  separa¬ 
tions  between  the  blue-  and  redshifted  masers  are  approximately 
860,  880,  and  890  AU,  respectively. 

In  order  to  better  characterize  the  net  expansion  of  the  ma¬ 
sers,  we  computed  the  separations  between  pairwise  combina¬ 
tions  of  components.  This  technique  has  previously  been  used  to 


TABLE  2 

FLO  Maser  Characteristics  Derived  from  the  Epoch  2  (2005  November  2)  VLBA  Data 


t>LSR  S„  as„  Relative  R.A.  or.a.  Relative  Decl.  <7Deci 

(kms-1)  (Jybeam-1)  (Jy  beam-1)  (mas)  (mas)  (mas)  (mas) 


Blueshifted  Masers 

-83.29 .  0.153  0.009  15.458  0.029  -10.516  0.045 

-81.57 .  1.323  0.010  0.000  0.017  0.000  0.028 

-81.41  .  0.096  0.010  -2.207  0.041  -0.456  0.063 

-80.73  .  0.216  0.010  8.174  0.032  0.173  0.051 


Redshifted  Masers 

-34.47 .  0.468  0.010  1.286  0.020  -112.898  0.033 

-33.54 .  0.293  0.009  2.518  0.020  -112.793  0.034 

-33.31  .  0.086  0.009  -6.432  0.048  -112.695  0.071 

-33.10 .  0.569  0.009  -2.658  0.020  -112.669  0.030 

-32.68  .  0.142  0.009  5.765  0.034  -110.732  0.051 
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TABLE  3 

H20  Maser  Characteristics  Derived  from  the  Epoch  3  (2006  February  24)  VLBA  Data 


nLSR  S„  crsLI  Relative  R.A.  crRA  Relative  Decl.  <TDeci 

(km  s-1)  (Jy  beam"1)  (Jy  beam-1)  (mas)  (mas)  (mas)  (mas) 


Blueshifted  Masers 

-83.29 .  0.174  0.013  15.580  0.029  -10.741  0.054 

-81.54 .  0.970  0.013  0.000  0.026  0.000  0.043 

-81.30 .  0.079  0.013  -2.358  0.063  -0.384  0.110 

-81.26 .  0.065  0.013  0.347  0.068  1.034  0.111 

-80.72  .  0.094  0.012  8.208  0.043  0.067  0.074 


Redshifted  Masers 

-34.44 .  0.228  0.015  1.315  0.039  -113.770  0.072 

-33.46 .  0.151  0.014  2.582  0.039  -113.587  0.070 

-33.35  .  0.061  0.012  -6.378  0.062  -113.447  0.106 

-33.04 .  0.164  0.012  -2.593  0.041  -113.582  0.075 

-32.71  .  0.123  0.012  5.717  0.047  -111.546  0.085 


characterize  OH  (Chapman  et  al.  1991;  Bloemhof  et  al.  1992; 
Kemball  1992),  H20  (Marvel  1996),  and  SiO  (Boboltz  et  al. 
1997;  Chen  et  al.  2006)  maser  motions  and  has  no  dependence  on 
the  alignment  of  the  maps  or  a  priori  assumptions  about  the  ve¬ 
locity  field.  The  procedure  involves  computing  the  angular  sep¬ 
aration  (A0)  between  two  features  at  one  epoch  (epoch  A)  and 
the  separation  between  the  corresponding  two  features  at  a  sec¬ 
ond  epoch  (epoch  B).  The  difference  between  the  two  values  of 
Ad  is  the  pairwise  separation  and  can  be  written  as 

A6B  -  A6a  =  \rt  -  rj\B-\ri  -  rj\A,  i=l,n;  j=i  +  l,n, 

(1) 

where  r,-  =  (x,-,  yt),  and  (xh  y,)  are  the  relative  offsets  in  right  as¬ 
cension  and  declination,  respectively.  The  procedure  is  repeated 


for  all  possible  pair  combinations,  and  the  separations  can  be 
plotted  as  a  histogram  (e.g..  Fig.  5).  The  inclusion  of  all  possible 
pair  combinations  often  results  in  a  bimodal  distribution  with 
one  of  the  peaks  biased  toward  zero.  This  is  because  some  of  the 
separations  involve  pairs  of  closely  spaced  maser  components  on 
the  same  side  of  the  distribution  (i.e.,  the  northern  or  southern 
regions  for  OH  12.8—0.9)  that  have  little  motion  relative  to  one 
another.  For  the  sake  of  clarity,  and  to  determine  representative 
values  for  the  angular  shifts  due  to  the  expansion,  we  have  in¬ 
cluded  only  those  pairs  separated  by  more  than  80  mas  in  the 
histograms  shown  in  Figure  5.  Figure  5  plots  the  pairwise  sep¬ 
arations  over  time  intervals  of  (a)  613,  ( b )  500,  and  (c)  1 14  days, 
respectively.  All  three  histograms  have  centroids  that  are  bi¬ 
ased  toward  positive  values  indicative  of  expansion.  The  mean 
(median)  differences  between  the  north  and  south  masers  are 


4  - 


0  - 


-12  - 


12 


8  4 

R.A.  Offset  (mas) 


^  1  1 

O  Epoch  1 

- 1 - ' - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 

OH12.8-0.9  South  O  Epoch  1 

□  Epoch  2 

-100 

□  Epoch  2 

A  Epoch  3 

A  Epoch  3 

-104 

- 

A 

w 

© 

f 

_ i _ 

cd 

s. 

• 

1  -108 

_ 

o 

•  •  •  • 

d 

- 

P 

■ 

▲ 

-112 

- 

o  m  ■  D 

*  i.  *  4 

-116 

_ 1 _ i _ 1 _ 

_ 1 _ i _ 1 _ l _ 1 _ i _ 1 _ l _ 1 _ 

4  0 

R.A.  Offset  (mas) 


-4 


-8 


Fig.  3. — Northern  blueshifted  masers  associated  with  OH  12.8—0.9  for  all 
three  epochs  (2004  June  21,  2005  November  2,  and  2006  February  24).  Masers 
are  color-coded  according  to  velocity,  and  point  size  is  proportional  to  the 
logarithm  of  the  maser  flux  density.  There  are  four  maser  features  common  to  all 
three  epochs,  including  the  reference  feature  at  (0,  0).  The  masers  show  little 
motion  relative  to  this  stationary  reference  feature. 


Fig.  4. —  Southern  redshifted  masers  associated  with  OH  12.8—0.9  for  all 
three  epochs  (2004  June  21,  2005  November  2,  and  2006  February  24),  dem¬ 
onstrating  the  motions  of  the  masers  over  time  relative  to  the  reference  feature 
in  Fig.  3.  Masers  are  again  color-coded  according  to  velocity,  and  point  size  is 
proportional  to  the  logarithm  of  the  maser  flux  density.  All  five  features  are  com¬ 
mon  to  all  three  epochs. 
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Separation  (mas) 


Fig.  5. — Histograms  showing  the  pairwise  spatial  separations  between  com¬ 
binations  of  northern  and  southern  maser  pairs  for  {a)  epochs  3  and  1 ,  ( b )  epochs  2 
and  1,  and  (c)  epochs  3  and  2.  Listed  in  each  panel  is  the  lapsed  time  between 
epochs  in  units  of  days. 

4.64  (4.69)  mas,  3.88  (3.89)  mas,  and  0.77  (0.80)  mas  for  613, 
500,  and  114  days,  respectively.  The  1  a  standard  deviations 
are  0.28, 0.21,  and  0.13  mas,  respectively.  The  equivalent  mean 
angular  velocities  are  2.76  ±0.16,  2.83  ±0.15,  and  2.47  ± 
0.41  mas  yr~ 1 .  To  within  the  errors,  the  three  velocities  are  con¬ 
sistent,  and  it  is  impossible  to  determine  whether  there  is  any 
acceleration  of  the  masers  in  the  outflow. 

If  8  kpc  is  again  assumed  to  be  the  distance  to  OH  12.8—0.9, 
the  linear  proper  motions  of  the  masers  can  be  computed  for  the 
angular  velocities  above.  The  resulting  average  linear  separation 
velocities  are  105  ±  6,  107  ±  6,  and  94  ±  16  km  s_1  for  the  613, 
500,  and  114  day  intervals,  respectively.  If  we  assume  that  the 
north  and  south  masers  are  moving  outward  at  equal  speeds  (i.e., 
at  half  the  computed  separation  velocity),  then  the  outflow  veloc¬ 
ity  in  the  plane  of  the  sky  would  be  ^53  km  s_1  for  the  613  day 
interval.  This  velocity  is  slightly  more  than  double  the  outflow 
velocity  detemiined  for  OH  12.8—0.9  in  Boboltz  &  Marvel  (2005), 
which  was  based  solely  on  the  radial  velocities.  In  the  following 
section  (§  3.2)  we  update  the  radial  velocities  of  the  masers  and 
combine  them  with  the  proper  motions  in  order  to  estimate  the 
full  3D  kinematics  of  the  masers  in  the  jet. 

3.2.  Three-dimensional  Maser  Outflow 

In  Boboltz  &  Marvel  (2005)  we  treated  the  radial  velocity  of 
OH  12.8—0.9  by  simply  determining  the  velocity  separation  of 
the  blue-  and  redshifted  maser  features  representing  the  peaks  in 
the  spectral  distribution.  This  was  an  improvement  over  previous 
single-dish  and  VLA  studies  that  contained  limited  spatial  infor¬ 
mation  and  likely  suffered  from  blending  of  features  in  the  spec¬ 
tral  domain.  In  light  of  the  two  new  epochs  of  VLBA  data,  we 
felt  that  the  more  rigorous  approach  of  computing  spectral  sepa¬ 
rations  (At>)  from  multiple  features  was  more  appropriate.  To  do 
this,  we  averaged  the  velocities  of  the  blue-  and  redshifted  ma¬ 
sers  identified  in  Tables  1-3  separately.  We  also  computed  the 
standard  deviations  for  each  set  of  components.  From  these  av¬ 
erages  we  computed  the  spectral  separation  Av  at  each  epoch. 


TABLE  4 

Velocity  Separation  of  Blue-  and  Redshifted  Water  Masers 


Obs.  Date  Av  a&v  Spectral  Res. 

(yr)  (km  s-1)  (kms'1)  (km  s-1)  Reference 

1987.452 .  38.8  1.5  0.16  1 

1988.516 .  37.7  1.6  0.16  1 

1990.132 .  42.0  2.1  0.16  1 

1990.249 .  41.2  2.3  0.16  1 

1991.332 .  41.5  2.1  0.16  1 

1992.049 .  42.0  2.0  0.16  1 

1992.975 .  42.3  2.3  0.16  1 

1994.184 .  40.9  2.1  0.16  1 

1995.186 .  42.4  2.0  0.16  1 

1996.675 .  43.3  2.3  0.16  1 

1999.088 .  44.6  2.0  0.16  1 

2004.473 .  48.7  1.6  0.22  2 

2005.838 .  48.3  1.3  0.22  2 

2006.151 .  48.2  1.3  0.22  2 


References. —  (1)  Engels  2002;  (2)  this  work. 

The  error  in  Av  is  simply  the  standard  deviations  added  in 
quadrature.  Using  this  method,  we  find  mean  values  for  At>  of 
48.7  ±  1.6,  48.3  ±  1.3,  and  48.2  ±1.3  km  s_1  for  epochs  1,  2, 
and  3,  respectively.  These  values  are  also  listed  in  Table  4. 

In  Boboltz  &  Marvel  (2005)  we  compared  the  spectral  separa¬ 
tion  of  the  blue-  and  redshifted  peaks  in  our  spectrum  with  peak 
separations  from  previous  single-dish  and  VLA  observations 
(i.e.,  Engels  et  al.  1986;  Gomez  et  al.  1994;  Engels  2002).  Here 
we  wish  to  compare  our  spectral  separations  computed  from  the 
component  averages  with  these  previous  results;  however,  only 
Engels  (2002)  tabulates  parameters  for  individual  spectral  fea¬ 
tures.  Engels  (2002)  categorizes  components  A-H  of  his  Table  4 
into  the  blueshifted  group  and  components  M-Q  of  his  Table  5 
into  the  redshifted  group,  with  the  remaining  components,  I-L, 
listed  as  intermediate.  From  this  information,  we  computed  av¬ 
erage  velocities  for  the  blue-  and  redshifted  spectral  features  and 
values  for  Ai>  in  a  manner  similar  to  that  performed  for  our  own 
VLBA  data.  The  resulting  values  of  Av  and  the  corresponding 
errors  for  the  Engels  (2002)  data  are  listed  in  Table  4.  The  only 
drawback  to  the  Engels  data  as  compared  to  our  VLBA  data  is 
the  fact  that  there  is  no  spatial  information;  thus,  the  degree  to 
which  spectral  blending  is  a  factor  is  unknown.  It  should  be 
noted  that  the  spectral  resolution  of  the  Engels  measurements  is 
slightly  better  than  our  VLBA  resolution  at  0.16  km  s  1 . 

Plotted  in  Figure  6  are  the  results  of  our  computation  of  Av 
from  the  blue-  and  redshifted  component  averages.  The  values 
from  the  1 1  epochs  of  Engels  (2002)  are  plotted  as  filled  trian¬ 
gles,  and  those  from  our  three  epochs  of  VLBA  data  are  plotted 
as  filled  circles.  The  data  clearly  show  an  increase  in  Av  as  a  func¬ 
tion  of  time.  To  these  data,  we  performed  a  linear  least-squares  fit 
weighted  by  the  squares  of  the  errors.  This  fit  is  plotted  as  the  line 
in  Figure  6,  with  a  slope  of  0.53  ±  0.04  km  s_1  yr-1  corre¬ 
sponding  to  the  relative  acceleration  of  the  blue-  and  redshifted 
masers  in  the  radial  direction.  This  is  slightly  lower  than  the 
value  of  0.68  ±  0.06  km  s  1  yr~*  reported  in  Boboltz  &  Marvel 
(2005).  With  the  new  procedure  of  averaging  over  multiple  spec¬ 
tral  features,  we  have  attempted  to  remove  the  error  involved  in 
computing  the  velocity  separation  from  single  peaks  on  the  blue- 
and  redshifted  sides  of  the  spectrum  that  are  likely  unrelated  over 
periods  greater  than  the  lifetime  of  the  components,  <3  yr  (Engels 
2002).  Engels  also  discusses  ~1  km  s_1  velocity  shifts  that  were 
observed  for  some  components  in  the  spectra  of  OH  12.8—0.9  and 
the  likelihood  of  these  shifts  being  caused  by  intensity  changes 
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Fig.  6. — Measured  velocity  separation  (Ar)  of  the  blue-  and  redshifted  water 
maser  features  as  a  function  of  time,  from  values  listed  in  Table  4.  The  line  rep¬ 
resents  a  weighted  linear  least-squares  fit  to  the  data  and  indicates  a  constant 
acceleration  of  0.53  km  s-1  yr-1. 

within  blended  features.  Such  shifts  could  affect  the  velocity 
separations  computed  from  single  spectral  peaks,  as  was  done  in 
Boboltz  &  Marvel  (2005). 

If  we  assume  an  average  over  our  three  VLBA  epochs, 
Av  «  48  km  s'1,  as  the  speed  at  which  the  masers  are  moving 
away  from  each  other  along  the  line  of  sight  to  the  observer,  then 
the  outflow  velocity  in  the  radial  direction  is  ~24  km  s_1.  The 
3D  outflow  velocity  of  the  masers  may  be  estimated  by  com¬ 
bining  the  radial  motion  above  with  the  motion  in  the  plane  of  the 
sky  determined  in  §  3.1.  We  find  the  3D  outflow  velocity  of  the 
masers  to  be  ~58  km  s'1,  with  an  inclination  angle  for  the  jet  of 


~24°  with  respect  to  the  plane  of  the  sky.  The  3D  outflow  ac¬ 
celeration,  assuming  this  inclination  angle  and  the  radial  accel¬ 
eration  above,  is  ^0.63  km  s'1  yr-1.  The  outflow  acceleration  in 
the  plane  of  the  sky  would  be  ^0.58  km  s' 1  yr' 1  under  this  as¬ 
sumption.  Thus,  over  the  613  days  between  epochs  1  and  3,  we 
could  expect  a  change  in  the  velocity  of  the  blueshifted  masers 
relative  to  the  redshifted  masers  of  1 .93  km  s'1  or  0.05  mas  yr'  , 
which  is  undetectable  in  our  present  observations. 

From  the  above  information  and  two  additional  assumptions, 
namely,  that  the  dynamical  center  of  the  outflow  is  the  midpoint 
along  the  axis  of  separation  between  the  two  maser  regions  and 
that  the  masers  have  zero  initial  velocity,  an  upper  limit  to  the  dy¬ 
namical  age  of  the  outflow  can  be  computed.  We  find  this  dy¬ 
namical  age  to  be  ^90  yr.  This  value  is  roughly  consistent  with 
the  age  determined  in  Boboltz  &  Marvel  (2005),  which  was  es¬ 
timated  with  no  knowledge  of  the  maser  motions  in  the  plane  of 
the  sky  or  the  inclination  angle  of  the  outflow.  If  one  assumes  that 
the  acceleration  has  and  will  remain  constant  at  the  above  value 
of  0.63  km  s'1  yr'1,  then  the  total  time  to  reach  an  outflow 
velocity  comparable  to  the  other  water-fountain  sources  (At;  ss 
150  km  s'1)  is  only  about  ~240  yr,  roughly  2.5  times  the  current 
dynamical  age  of  OH  12.8—0.9. 

In  order  to  compare  the  characteristics  of  OH  12.8—0.9  to 
those  of  the  other  water-fountain  sources,  we  have  summarized 
the  relevant  properties  of  all  four  sources  in  Table  5.  From  the 
table  we  see  that  the  assumed  distance  to  OH  12.8—0.9  is  inter¬ 
mediate  among  the  four  sources.  The  peak  flux  density  is  shown 
as  a  range  measured  over  multiple  epochs  by  Engels  (2002)  for 
OH  12.8—0.9  and  by  Likkel  et  al.  (1992)  for  the  three  other 
sources.  The  peak  flux  range  for  OH  12.8—0.9  is  similar  to  that 
of  the  other  distant  source  IRAS  191 34+2 131  and  weaker  than 
those  of  W43Aand  IRAS  16342—3814.  The  angular  extent  of 
the  outflow  on  the  plane  of  the  sky  is  also  similar  to  that  of 
IRAS  19134+2131.  The  linear  extent  of  OH  12.8— 0.9,  however, 
is  much  less  than  those  of  any  of  the  other  water-fountain  sources. 
An  increase  in  the  assumed  distance  to  the  source  would  serve  to 
bring  this  value  in  line  with  the  other  sources.  For  the  collimation 
of  the  OH  12.8—0.9  jet,  we  have  used  an  average  of  the  northern 
and  southern  arcs  and  find  that  the  collimation  is  slightly  larger 


TABLE  5 

Summary  of  the  H+)  Maser  Characteristics  for  the  Water-Fountain  Sources 


Characteristic  OH  12.8-0.9  W43A  IRAS  19134+2131  IRAS  16342-3814 

Distance  (kpc) .  8  2.6  16  2 

«W*(Jy) .  t-4  13.0  11.5-56.6  1.2-6.0  16.0-67.7 

Angular  extent13  (mas) .  110  920  135  3000 

Linear  extent13  (AU) .  880  2400  2180  6000 

Outflow  collimation  (deg) .  15  5  10  6 

Inclination  anglec  (deg) .  24  39  25  40d 

Outflow  velocity: 

vtan  (km  s-1) .  53  110  120 

iijadfkms"1) .  24  90  65  130 

vtot  (km  s-1) .  58  145  130 

Outflow  acceleration: 

tfrad  (km  s-1  yr-1) .  0.53 

Dynamical  age  (yr) .  90  50  50  150d 

References .  1, 2,  3  4,  5,  6  4,  7  4,  8,  9,  10 


a  Range  over  multiple  epochs  of  single-dish  observations. 
b  Measured  in  the  plane  of  the  sky. 
c  Relative  to  the  plane  of  the  sky. 
d  From  OH  maser  observations. 

References. — (1)  Engels  2002;  (2)  Boboltz  &  Marvel  2005;  (3)  this  work;  (4)  Likkel  et  al.  1992;  (5)  Imai  et  al.  2002;  (6)  Imai 
et  al.  2005;  (7)  Imai  et  al.  2004;  (8)  Sahai  et  al.  1999;  (9)  Morris  et  al.  2003;  (10)  Claussen  et  al.  2004. 
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than  that  for  the  other  three  water-fountain  sources.  If  we  use  only 
the  southern  arc,  or  if  we  disregard  the  one  northern  feature  at 
(15.0,  —10.0),  then  the  collimation  is  ~  1  O  ',  which  is  consistent 
with  the  other  water-fountain  sources.  The  estimated  inclination 
angle  of  the  outflow  is  also  comparable  with  values  estimated  for 
the  other  three  objects. 

Aside  from  the  difference  in  linear  extent,  the  primary  differ¬ 
ence  between  OH  12.8—0.9  and  the  other  water-fountain  sources 
is  in  the  kinematics.  Both  the  radial  (t>rad)  and  tangential  (t>tan) 
outflow  velocities  are  less  than  half  the  values  for  the  other  sources. 
It  is  possible  that  the  tangential  velocity  could  be  closer  to  that 
of  the  other  sources  if  the  distance  to  the  source  turns  out  to  be 
greater  than  8  kpc.  An  increased  distance  would  not  increase  the 
radial  velocity  estimate,  but  would  instead  result  in  a  decreased 
estimate  for  the  inclination  angle  of  the  jet.  Another  major  differ¬ 
ence  between  OH  12.8—0.9  and  the  three  other  water-fountain 
sources  is  that  OH  12.8—0.9  appears  to  be  the  only  object  with  a 
measurable  acceleration  (arad)  of  the  masers  along  the  line  of 
sight.  Likkel  et  al.  (1992)  observed  the  H20  maser  spectra  of  the 
three  other  water- fountain  sources  over  the  course  of  a  few  years 
and  found  no  evidence  for  systematic  velocity  drifts  in  spectral 
features  that  would  suggest  acceleration  of  the  masers.  The  ac¬ 
celeration  in  combination  with  the  smaller  extent  of  the  masers 
toward  OH  12.8—0.9  might  suggest  that  it  is  younger  than  the 
other  three  objects.  This  is  not,  however,  confirmed  by  the  com¬ 
puted  dynamical  age  of  the  jet,  which  is  intermediate  among  the 
four  sources  at  90  yr. 

One  characteristic  that  is  not  shown  in  the  tables  is  the  mor¬ 
phology  of  the  masers  relative  to  the  central  star  and  the  jet  di¬ 
rection.  From  the  images  of  OH  12.8—0.9  presented  here,  we 
find  that  the  masers  are  arranged  in  point-symmetric  fashion  along 
arcs  perpendicular  to  the  direction  of  the  outflow,  suggesting  that 
the  outflow  is  impacting  a  denser  surrounding  medium.  This  bow 
shock  structure  is  also  observed  for  the  H20  masers  toward  IRAS 
16342—3814,  at  least  for  the  blueshifted  side  of  the  jet  (Claussen 
et  al.  2004).  These  two  sources  are  in  contrast  to  W43A,  where 
the  masers  are  arranged  along  the  direction  of  the  outflow.  The 
structure  of  the  W43  A  masers  is  well  represented  by  a  precessing 
jet  model  (Imai  et  al.  2002, 2005).  The  remaining  water-fountain 
source,  IRAS  19134+2131,  shows  the  point-symmetric  mor¬ 
phology  of  the  three  other  sources,  but  the  relationship  between 
the  masers  and  the  outflow  is  less  clear.  From  the  two  epochs 
of  VLBI  observations  reported  in  Imai  et  al.  (2004),  the  masers 
on  the  eastern  redshifted  side  appear  to  be  roughly  aligned  with 
the  line  connecting  them  to  the  western  blueshifted  masers.  The 
motions  of  the  masers,  however,  do  not  appear  to  be  along  this 
line,  but  rather  perpendicular  to  this  in  a  primarily  north-south 
direction. 

Similar  to  the  W43  A  outflow,  precession  is  indicated  for  IRAS 
16342—3814,  but  not  as  a  result  of  the  available  H20  maser  data. 
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Instead,  it  was  the  Keck  adaptive  optics  images  in  the  near- 
infrared  (Sahai  et  al.  2005)  that  showed  a  corkscrew  structure  in¬ 
scribed  on  the  walls  of  the  observed  lobes.  The  combined  radio/ 
infrared  observations  demonstrate  that  both  bow  shock  morphol¬ 
ogies  and  precession  can  take  place  in  the  same  object.  Since 
no  such  high-resolution  optical /infrared  imaging  is  available  for 
OH  12.8—0.9,  it  is  unknown  whether  precession  occurs  in  the  ob¬ 
ject;  however,  it  is  clearly  not  suggested  by  the  structure  or  kine¬ 
matics  of  the  H20  masers. 

4.  CONCLUSIONS 

Using  the  VLBA,  we  have  shown  that,  like  the  spatial  mor¬ 
phology,  the  kinematics  of  the  H20  masers  toward  OH  12.8—0.9 
indicate  that  the  object  is  a  member  of  the  rare  “water-fountain” 
class  of  sources.  We  find  that  the  masers  continue  to  be  located 
in  two  distinct  regions  at  the  ends  of  a  bipolar  jetlike  structure 
having  a  north-south  orientation  on  the  sky.  The  masers  in  the 
two  regions  are  distributed  along  arclike  structures  ~12-20  mas 
across  that  are  oriented  roughly  perpendicular  to  the  separation 
axis.  The  angular  separation  of  the  two  regions  is  roughly  110  mas, 
with  the  blueshifted  masers  located  to  the  north  and  the  redshifted 
masers  to  the  south.  The  arclike  arrangements  of  the  masers  are 
suggestive  of  bow  shocks  formed  by  a  collimated  axisymmetric 
wind  impinging  on  the  ambient  medium  surrounding  the  star. 

The  two  additional  epochs  of  VLBA  observations  beyond  our 
first  epoch  in  2004  have  allowed  us  to  track  the  projected  spatial 
motions  of  individual  features  over  the  course  of  ~600  days.  We 
find  that  the  two  maser  regions  maintain  their  arclike  appearance 
as  they  expand  away  from  each  other  along  the  axis  of  separa¬ 
tion.  The  observed  proper  motion  of  the  southern  masers  relative 
to  the  northern  masers  is  2.76  mas  yr-1  (^105  km  s~!  at  the 
assumed  distance  of  8  kpc).  The  masers  maintain  collimation 
during  this  expansion,  with  little  motion  perpendicular  to  the 
north-south  axis.  Combining  the  proper  motions  with  the  mea¬ 
sured  radial  velocities  yields  a  three-dimensional  (3D)  expansion 
velocity  of  58  km  s~[  relative  to  the  dynamical  center  of  the  dis¬ 
tribution.  The  inclination  angle  of  the  jet  is  estimated  to  be  ~24°. 
A  simple  linear  fit  to  the  historical  radial  velocity  data  combined 
with  our  own  measurements  yields  an  estimate  for  the  3D  accel¬ 
eration  of  ^0.63  km  s~ 1  yr~ 1  and  a  corresponding  dynamical  age 
for  the  outflow  of  ~90  yr.  Long-term  VLBA  monitoring  of  the 
H20  masers  should  enable  the  determination  of  the  tangential 
acceleration  of  the  masers  and  the  true  3D  acceleration  of  the 
outflow  associated  with  OH  12.8—0.9. 
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